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Abstract 

In this series of lectures we review observational evidence for, and theoretical investiga- 
tions into, cosmic acceleration and dark energy. The notes are in four sections. First I 
review the basic cosmological formalism to describe the expansion history of the universe 
and how distance measures are defined. The second section covers the evidence for cos- 
mic acceleration from cosmic distance measurements. Section 3 discusses the theoretical 
avenues being considered to explain the cosmological observations and section 4 discusses 
how the growth of inhomogeneities and large scale structure observations might help us 
pin down the theoretical origin of cosmic acceleration. 

1 Describing spacetime: distances and General Rel- 
ativity 

Here we review the basics of General Relativity that allow us to describe the expansion 
history of a homogeneous and isotropic Universe - a reasonable description of the Universe 
on the largest cosmic scales. 

1.1 Units 

First a note on units. Throughout, we will often set c = ks = h = 1. This is just a change 
of units making time, distance and temperature consistent with units of energy, 

h X time~^ •<->■ he X distance'^ fcs x temp energy (1) 

It is often useful to use astrophysical, as well as SI, units of length. 

• The Astronomical unit is a good unit for solar system scale distances: 

Astronomical unit(AU) = Earth — Sun distance (2) 

= 1.49 X 10"m (3) 
Solar system ~ AOAU (4) 

• Trigonometric parallax (how the apparent position of objects moves as the Earth 



orbits the Sun) defines the parsec unit of length 



Parsec(pc) = Distance to object with parallax of 1 arcsec 

lAU 

Ipc = 

larcsec 

= 206,265AC/ 
= 3.09 X lO^^m 

The pc is useful unit to describe astrophysical distances 
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1.2 Cosmic dynamics 



On scales above a few hundred Mpc the universe is well-approximated as homogeneous 
(no preferred location) and isotropic (no preferred direction). Being homogeneous and 
isotropic implies 

• the universe can have one of three types of geometry: negatively curved, positively 
curved or flat, with a radius of curvature, that is the same at every point in 
space. 

• the universe's expansion rate has to be the same at every place, i.e. it has only 
time, but no spatial, dependence, described by a scale factor a{t). We will assume 
a convention where a = 1 today. 

The relationship between distances and time intervals in a homogenous and isotropic 
universe can be conveniently described by the Friedmann Robertson Walker (FRW) 
metric. This is commonly written in one of two, equivalent, forms: 

Using spherical polar coordinates, r is the comoving radial distance, (Kl — dO^ + 
sin^ OdcjP' and physical time 



ds'^ = -de + a{tf [dr'^ + Sl{r, R)dQ] 
or with a change of variable x = S^ir, R), K = n/R^, 

dx^ 



ds' 



-dr + a(ty 



1-Kx^ 

where k, K = k/R^ and S^iriR) depend upon the geometry 



(13) 



(14) 
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Flat (K=0) 

Paths diverge at constant angle 
Angles in triangle =180° 



Positive curvature (K>0) 

Paths diverge 

Angles in a triangle >180° 



Negative curved (K<0) 
increasingly converge 
Angles in a triangle <180° 




Figure 1: The three types of curvature consistent with a homogeneous and isotropic space. 



We can also define conformal time dr = dt/a{t) so that the FRW metric can be written 



Finally we can write the space time element, ds"^, in a coordinate independent way, 
in terms of a metric, g^^, for a generic 4D space-time coordinate system x'^ where fi = 
0, 1, 2, 3. fi = is reserved for the time coordinate and /x = 1, 2, 3 are for the 3 spatial 
coordinates. If we are just referring to the spatial coordinates then one typically uses 
Roman letters for the indices, e.g. gij where i,j = 1, 2, 3. 

In Einstein's theory of General Relativity, the effect of gravity is to distort the cur- 
vature of space in the presence of matter. The relationship between the curvature and 
matter is given by Einstein's equation 



The left hand side of the equation describes the curvature of space: is the Einstein 
tensor, R^^, is the Ricci tensor and R = g^^R^i, is the Ricci scalar, both of which are 
functions of the metric g^u and its derivatives with respect to the space-time coordinates. 
The RHS is written in terms of T^j,, the energy-momentum tensor, describing the 
properties of matter in the universe. 

For a perfect fluid, with energy density p, and isotropic pressure, P, 




(15) 



Gfj,,y = Rfj,u — -g^iuR — SvrGT^jy. 



(16) 



T^''' = {p + P)U^'U'' + Pg- 



(17) 
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where If^ is the four velocity of the fluid. In the fluid's rest frame [/'^=diag(l, 0,0,0) and 



f P P P\ 

T,. = diag^p,^,^,^J (18) 
= diag (-p,P,P,P) (19) 
The n = u = component of Einstein's equation gives us the Priedmann equation 

SnG K , , 

= ^ = — '2°' 

where dots denote derivatives with respect to t. 

The Friedmann equation tells us that a universe containing matter has to be dynam- 
ically evolving, either expanding or contracting, it cannot be static. Commonly one can 
also find Newton's constant replaced by the reduced Planck mass = I/SttG (in units 
where h = c = 1). 

The combination of the fi = u = j and /i = z/ = equations gives us the acceleration 
equation 

^ = -^^(P + 3P). (21) 
a 6 



For "normal" matter the energy density and pressure are both positive for which (21) 
predicts decelerating cosmic expansion. 

Energy-momentum conservation allows us to see how the energy density of matter 
evolves as the universe expands, 

rpfl _ rjifi I -p^ rjia _ "pO ™ _ r| ( cycy\ 

v,^ ~ afi-'- ij, i/^-i- a "J- \^^/ 

Here = dX/dx^, a partial derivative with respect to the coordinate x'^, X.^ = V^X 
denotes the covariant derivative, that preserves parallel transport in a curved or expanding 
space and F^^ is the Christoffel symbol relating the two derivatives. In flat non-expanding 
space, all F^^ = 0, and the two derivatives are identical. The covariant derivative is such 
that o^r=0. 

For u = we get the fluid or "continuity" equation, 

p + 3H{p + P) = (23) 

Two out of these three equations (the Friedmann, acceleration and fluid equations) 
are independent, the third being able to be obtained from the other two. 



1.3 The Cosmological Constant 

In classical, non gravitational physics, the normalization of energy is arbitrary, since the 
dynamics of particle moving under some potential, V{x), is the same as that moving 
through, V{x) + Vq. In GR however the actual value of the energy density matters, 
i.e. the value of Vq has a role in the dynamics. This opens up the freedom to include 
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a cosmological constant, A, a component with a constant energy density, into the 
dynamical equations. Here we discuss the effect of including A to the dynamical equations, 



and reserve discussions of its theoretical origins to section 3.1 



If the effect of adding in a constant energy density is to be insensitive to the choice of 
coordinates then it has to satisfy 



rptJ.{A) 







(24) 



from our discussion in the previous section about the effect of the covariant derivative on 



the metric, this implies that T'^j^ oc 5^'^^. This means that from (19), T^^ = —p^Q^u and 
Pa 



-p\. Note Pa = — pA also falls out of the fluid equation (23), requiring pa = 



Including A, as an interpretation as a new type of matter, modifies the RHS Einstein 
equation 



(25) 



Including A does not destroy the general covariance (coordinate independence) of 
Einstein's equations. For this reason Einstein was able to include A to the LHS of his 
equations, interpreting it as a mathematical degree of freedom that can be added. 



G^u + ^Oi^iy = StiGT^ 



(26) 



Both of these are mathematically and dynamically equivalent. 

Adding in a cosmological modifies the Friedmann and acceleration equations 



SvrG 
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■^(P + 3P) + - 



(27) 
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This allows the universe's expansion to be accelerating if A dominates. 
1.4 The equation of state parameter, w 

For a matter species i we can describe the relationship between its pressure and energy 
density using an equation of state parameter, Wi, 



pi 



We can write the fluid equation, (23), in terms of this 



which when integrated yields 

Pija) 



p + 3H{l + w)p = 



exp 



-3 / (1 + Wi{a))d\na 

J a 
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(29) 



(30) 



(31) 



where p° is the energy density today, when a = 1. 

For normal matter, the pressure of a species is related to the kinetic energy oc m < 
V >^, while the energy density ~ mc^, hence for non-relativistic matter, such as baryons 
and cold dark matter (CDM), one would expect the equation of state to be negligible, 
while for relativistic species, like photons, < v"^ >~ 

We can describe these, and A, very well using a constant equation of state 



Wnon-rel ~ (32) 
Wrel = 2 (33) 
WA = -1. (34) 



For a constant w (23) yields 



a-'^'+'"^> (35) 



The universe is made up of a variety of different matter species, and to consider the 
dynamics from their combined effect we can define an effective equation of state, 

Pt t 

wtot = — (36) 

Ptot 



which from (21) will give acceleration if Wtot < 1/3. 



1.5 Dark energy 

The cosmological constant can generate acceleration with w = —1, however, as we will 
discuss in detail in subsequent lectures, it is not alone in this. The family of theoretical 
types of matter that can give rise to acceleration are given the collective name dark 
energy. They are typically described by an equation of state parameter w that is either 
constant, or time evolving, but that is sufficiently negative to give Wtot < ~l/3 at late 
times. 

1.6 The critical density, Pcrit^ and fractional energy density, Qi 
The critical density, Pcrit{o)i is the total energy density of matter required to give the 



Hubble parameter H{a) in a fiat universe, i.e. = in (20) 



_ 3H\a) 

PcriM) = g ^ (37) 

The fractional energy density in species f2j(a),is the fraction of the critical density 
in that species 

a(a) = (38) 

The fractional energy today for a species will be written 
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1.7 The deceleration parameter, q 



To complement the Hubble parameter, H{a) = a/ a we can define an acceleration pa- 
rameter, q{a), to describe the second order change in the expansion history, 

,(a) . - J,. (39) 
If GR is assumed to hold then the Friedmann and acceleration equations give 

q{a)^lj2^,{a)[l + 3w,{a)], (40) 



1.8 Redshift, z 



The redshift, z, measures the stretch in the wavelength of light due to the expansion 

of space. Consider two photons emitted at time and tg + dtg, and received at tr and 
tr + dteT. Since they have traversed the same comoving distance 



/•te+dte 
Jte 



dt' 

ait') 


ptr^dtr 
J te-^-dte 


dt' 

ait') 


(41) 


dt' 

a{t') 




dt' 

ait') 


(42) 


dte 

a(te) 


dtr 

a{tr) 




(43) 


X{tr) 


a{tr) 
a{te) 


1 
a 


(44) 



Define the redshift 

X{te) a 

Note this cosmological redshift is caused by the expansion of space, not due to the motion 
of the galaxies that would exist in the absence of expansion (this is called peculieir 
motion). 

(1 + Ztot) = (1 + Zcosmo)ii + Zpec) (46) 

Typical peculiar velocities are of order a few 100 km/s so that, given Hq, peculiar velocities 
are ~ 10% of cosmological expansion rate at distances ~ few tens of Mpc. 



1.9 Comoving distance, x 

The comoving distance, x is the distance between two objects instantaneously today. 
It is impossible to actually measure this, because of finite speed of light, but nonetheless 
it is a useful conceptual distance to calculate since other, observable, distances can be 
simply related to it. 
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Consider a light ray purely in the radial direction dx = dr, d9 = dcj) = 0, then 



r 



X= dr (47) 



We can relate the distance to the path of a photon, and hence time, using the metric (13). 
For a photon, ds^ = 0. If it is emitted by the object at time t and observed today, at time 

a{t)dx = cdt (48) 
df 

Xit) = / ^ (49) 
We can alternatively write this in terms of Hubble factor 

The comoving distance can then be calculated for a given cosmology using the Friedmann 
equation 

The comoving horizon is the maximum distance, on the comoving grid that infor- 
mation can have been carried. It is the comoving distance light could have travelled since 
the start of the universe. 

da' , , 

(52) 



1.10 Angular diameter distance, dA 

The angular diameter distance, rf^ is the distance to an object inferred from comparing 
the angular size of an object to its (assumed known) length. 

We define the angular diameter distance using the naive Euclidean (flat), non-expanding 
geometry. If the object of length L subtends an angle 

dA = - (53) 
a 

We have to take into account that space has expanded between emission and observation. 
If light was emited at scale factor a{t), then the comoving size today is I /a. The angle, 
a, is unaffected by expansion. The comoving distance to the object is, x{(^)j comoving 
coordinates (i.e. length scales today) is then 

L , , 

a = — - 54 
ax{a) 
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hence 

dA = ax{a)^^ (55) 
i + z 

If we live in a curved space then the angular size of the object is modified by the curvature, 
through its effect on the path of the light rays from the object. Consider the object length 
L spans coordinates (r, 6, 0) and (r, 6 + 66, 0), then the physical length of the object in a 
generally curved space, from the RW metric is: 

L = a{t)S^{r,Ro)6e (56) 

In a curved space therefore 

dA = aS^{x,Ro) (57) 



1.11 Luminosity distance, di 

The luminosity distance, di is inferred from comparing the inferred brightness (or 
"flux" measured in power per unit area) of an object to its (assumed known) luminosity 
(in units of power). 

We define the luminosity distance again using a (naive) Euclidean, non-expanding 
geometry. If the object of luminosity L is observed to have flux F, then using the inverse 
square law. 



Cosmic expansion has two primary effects on the luminosity distance 1) the light is 
redshifted, the wavelength increases ~ 1/a and therefore the energy decreases by factor 
a, 2) the number of photons arriving at the detector is also reduced by a factor a because 
of the stretching of space, therefore Lof,s — a^Lemit- The comoving distance to the object 
is again xi'^)^ hence 



^TTxiay 



(59) 



= ^ = (1 + z)x{z) (60) 
a 

Consider photons emitted from a surface area A. Then if the universe's geometry is 
curved the effective observed area is AS'^{r, Rq) i.e. photons are bent outwards, increasing 
the effective area, decreasing the observed flux, in closed geoemetries and the opposite for 
open geometries: 

d, = ^^^^^ (61) 

A key difference from our naive Euclidean perspective is that different methods of 
measuring distances yield different results in an exanding universe. 

dL{z)^{l + zfdA{z) (62) 

While this effect becomes neghgible at low redshifts it is important at cosmologically 
relevant scales. 
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2 Evidence for cosmic acceleration from distance mea- 
surements 

2.1 Super novae 

Type la supernovae (SNIa) gave the first evidence for cosmic acceleration, and have 
arguably provided the most effective constraints on the dark energy equation of state 
parameter. SNIa are believed to be carbon-oxygen white dwarf stars in a binary system 
that are accreting matter from a companion. When the white dwarf's mass approaches 
the Chandrasekhar mass, ~ 1.4 solar masses, it undergoes a thermonuclear explosion. 

An object is a standard candle if its luminosity can be inferred from either an 
understanding of its physical properties, or some other observable. Astronomers typically 
refer to magnitudes rather than fluxes and luminosities when presenting observational 
data. An object's apparent magnitude, m, is related logarithmically to its flux. For 
two objects, their apparent magnitudes and fluxes are related by 

F2 

mi -1712 = 2.5 Ig— (63) 

with the zero-point m = set by a reference star, Vega. The Sun's apparent magnitude 
is rrisun = —26.7. An object's absolute magnitude is related to the absolute luminosity 
of the object. It is defined as its apparent magnitude if it were a fixed distance, lOpc, 
away from us. Using the definition of the luminosity distance 

m-M = 51g— . (64) 
lOpc 

The Sun's absolute magnitude is 4.83. 

The lightcurve for a SNIa explosion (the apparent magnitude m, versus time) is 
predominantly determined by radioactive decays of two species. The decay of ^^Ni dom- 
inates in the first few days/weeks, and determines the peak luminosity of the event. The 
lightcurve decay, in the latter weeks, is determined by the decay of ^^Co. If the white 
dwarf combustion is complete then it is expected that roughly ~ 0.6 solar masses of 
^^Ni is produced. As such, one might expect all SNIa events to have comparable lu- 
minosities and for them to be strong standard candle candidates. Since the typical 
peak luminosity of a SNIa is a few billion times that of our Sun [Msn ~ —19.3, im- 
plying LsN/Lsun = 10('^^^"-^sn)/5 ^ iq9 _ iqw^ ^ thesc events can be as bright as their 
host galaxy, providing an effective distance measure to large cosmic distances. In reality 
there is an intrinsic variation in the peak luminosities of SNIa which could threaten their 
usefulness as a standard candle. However an empirical relation between the SNIa peak 
brightness and the rate at which the lightcurve decays after the peak has been found, 
the Phillips relation [1], in which intrinsically brighter supernovae decay more slowly. 
This relation allows the luminosity to be inferred from fitting the lightcurve, reducing the 
intrinsic spead. 



From (58), the comparison of the observed apparent magnitude, m, to the presumed 



absolute magnitude yields an estimate of d^. Consider a Taylor expansion of the expansion 
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rate close to today (time t = to), 



a{t) = 1 + Ho{t - to) - ^(t - to)' + ... (65) 



where go is the value of the deceleration parameter (39) today. By writing a{t) in terms 



of [1 + z) ^ and inverting the expression to find t{z), one can obtain a Taylor expansion 



ofdL (jeoj), 

1 



HodLiz) = z + ^{1 - qo)z' + ... (66) 



From (40), go = — ^^a ^oi a ACDM universe and a measurement of go would lead 

to a degeneracy in the {Qm,^A} parameter space described by — = constant. 

One wouldn't be able to constrain Qm and f^A separately. Of course the Taylor expansion 
breaks down at higher z so that the constraint is not exactly of the form we see here, 
however a degeneracy still persists, and one needs other, complementary observations to 
break it. 

As an example, let's consider 2 supernovae: 

1992P : z = 0.026, m = 16.08 (67) 
1997ap: z = 0.83, m = 24.32. (68) 

For z ^ 1 then ^^(z) ~ z/Ho, and assuming a value for Ho = 72km/s/Mpc~^, one 
finds the absolute magnitude of the low z SNIa to be M = —19.09. Assuming this also 
holds for the higher supernova gives Hodilz = 0.83) ~ 1.16. One can then compare 
this to the theoretical estimate for Hodi for a pure matter universe Qm = 1, for which 
HodLiz = 0.83) ^ 0.95, and for a ACDM universe, fi^ = 0.3, ^a = 0.7, for which 
HodLiz = 0.83) ^ 1.23. 

Observational results for luminosity distances are often given using the distance mod- 
ulus, fi, 

PL = m-M = 5lg{HodL) (69) 

with 

M = M -5\g{HoMpc-^) + 25. (70) 

where A4 is calibrated using low-redshift supernovae. 

In reality, the measured apparent magnitude can be affected by numerous systematic 
errors including extinction and redenning from dust in the galactic host. The magnitude 
is measured over a finite frequency range, or passband. The distribution of the magnitude 
over the different passbands is referred to as the supernovae's color. Redenning can cause 
photons to be shifted from one passband to another which can lead to over or under- 
estimation of the apparent magnitude if a limited number of passbands are used. If 
the intrinsic colors of the SNIa or the extinction law were precisely known, then the 
extinction can be eliminated from the distance modulus by measurements of multiple 
passbands. Realistically however, intrinsic variations in galactic dust properties and SN 
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colors twinned with errors in the measurements of the apparent magnitude (photometric 
errors) will still lead to uncertainties in the distance modulus. A way to mitigate this is 
to measure the supernovae in the near-infrared, where the effects from dust extinction are 
not as significant. 

The accuracy of the distance modulus is also dependent on the low-redshift sample of 
supernovae used to "anchor" the Hubble diagram, and calculate Ai. Until recently this 
sample has been rather small and heterogenous and has correlated large-scale, correlated 
peculiar velocities. The peculiar velocities obscure extraction of the cosmological velocity 
and lead to errors in the distance estimates. Currently a number of surveys are under- 
way to increase the low redshift sample, including the Harvard Center for Astrophysics 
supernovae project [2], the Carnegie Supernova project, the Nearby Supernova Factory 
and the Sloan Digital Sky Survey II supernova survey. 

As supernovae at higher and higher redshift are observed the issue of evolution in 
supernovae luminosities arises. There is an observed correlation between peak bright- 
ness and host-galaxy type, and with host-galaxy properties (such as star formation rate 
and metalicity, and redshift). This evolution sensitivity can be reduced by using Phillips- 
corrected lightcurves are less sensitive to galaxy-host environment and finding low-redshift 
analogs of higher-redshift galactic environments to calibrate supernovae behavior. Ulti- 
mately though a better understanding of supernovae physics will be required to reduce 
sensitivity to supernovae evolution and the effect of environment. 

2.2 Cosmic Microwave Background (CMB) 

The Cosmic Microwave Background (CMB) are relic photons released at the time of 
recombination, when the ionized plasma of free electrons and protons present in the 
early universe cools enough for neutral hydrogen to form. While photons readily scatter of 
free electrons through Thomson scattering they interact negligibly with neutral hydrogen, 
so post-recombination the photons free-stream, with minimal interaction until they are 
detected today. 

The CMB photons are remarkably homogeneous in temperature across the whole sky. 
However small fluctuations in the temperature are present at the level of 10~^, or a few fiK 
in magnitude. The utility of the CMB in placing constraints on dark energy comes from 
considering the characteristic scale over which these fluctuations should be correlated. 

The CMB provides a complementary measure of the expansion history to the super- 
anovae. While supernovae are standard candles, the CMB is said to be a standard 
ruler. The angular diameter distance is inferred by comparing the apparent angular 
size of correlations in the CMB temperature to their expected intrinsic size, estimated 
from understanding the propogation speed of photons in the ionized plasma, prior to 
recombination. 

Recombination takes place when the photons temperature is T ~ 30007^'. As the 
universe expands the photons redshift and cool, and have a measured temperature today of 
T = 2.726i^. This implies that recombination occured when the universe was roughly 1100 
times smaller than it is today, Zrec ~ 1100. For the standard cosmological composition, 
this was when the universe was ~ 400, 000 years old. 

The speed of photons in the plasma prior to recombination, the sound speed Cg, is 
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roughly given by Cs = c/a/3. At recombation the sound horizon, r^, describes the max- 
imum distance photons could have travelled, and hence a typical correlation length for 
any inhomogeneities (temperature fluctuations) that could be generated through causal 
processes. The comoving scale of the sound horizon is given by 

which is ~ 145Mpc for a standard cosmological composition. Since recombination occurs 
early in the universe's history one finds that curvature and, in most scenarios, dark energy 
have negligible effects on the sound horizon. The apparent angular size of the correlations. 



6, is then determined to the angular diamater distance to recombination, following (53) 



The universe's curvature/ geometry could have a significant effect on the angular size 
these correlations subtend, with an open geometry decreasing 6*, and a closed geometry 
increasing it. The main constraint on the {VLm-, parameter space coming from the 
CMB is therefore a constraint on Vttot = + ^a, complementary to the measured 
combination from the supernovae. Since late time acceleration increases the angular 
diameter distance, adding dark energy or A has a geometric degeneracy with curvature. 
The degeneracy can be broken however by including an additional constraint on the value 
of ffo. 

The CMB contains more information than just the angular diameter distance that is 
sensitive to the properties of dark energy. This is from the Integrated Sachs Wolfe (ISW) 
effect related to the growth of large scale structure. We will discuss this later in section 



2.3 Baryon Acoustic Oscillations (BAO) 

Prior to recombination the photons and baryons (mostly protons) were tightly coupled to 
one another, through Thomson scattering, and the electrostatic attraction of the electrons 
and protons. Fluctuations in the baryon density were able to propogate at comparable 
speeds to the photons, i.e. with Cg ~ c/a/3. When recombination occured the CMB 
photons "decoupled" from matter. Since the baryon mass is much larger than their 
thermal energy the particles are highly non-relativistic and their sound speed dropped 
effectively to zero. After decoupling, a correlation with characteristic size remains 
imprinted in the distribution of baryons. Figure |2] gives a pictorial summary of how the 
characteristic scale is imprinted in the baryon and dark matter matter distribution we 
observe today. 

Fluctuations in the baryon density seed the galaxies and clusters of galaxes, large 
scale structure, and the galaxy correlation function that we observe today contains 
a preferred comoving scale ~ r^. When Fourier transformed, the correlation gives rise 
ripples in the power spectrum of the baryon density fluctuations. This correlation is 
therefore known as baryon acoustic oscillations. 
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Figure 2: Figures i-vi show the evolution of a delta-function spike in the baryon (labeled 
"gas") density over time along with analogous evolution for dark matter, photon and 
neutrinos. Figure i: An initial delta-function spike, early in the universe's history, well 
before recombination occurs. Figures ii and iii: the photon and gas over-densities are 
tightly coupled up to recombination, with both propogating at speed c^. At decoupling 
the correlation length (here depicted by the comoving radius of the overdensity "ripple" ) 
is ~ IbOMpc. Figure iv: following decoupling, the acoustic oscillation correlations in the 
baryon density distribution do not propogate further and the radius of the correlation 
remains fixed while the photons continue to propogate at the speed of light. Figure 
v: the dark matter and baryons now fall into the other's gravitational well. Figure vi: 
today we see correlated density distributions of baryons and dark matter that include 
the correlation scale imprinted from the baryon photon coupling prior to recombination. 
[Credit: Daniel Eisenstein and the Sloan Digital Sky Survey collaboration.] 
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Figure 3: The galaxy correlation function, ^, from the Sloan Digital Sky Survey |3] clearly 
showing the baryon acoustic correlation at 100h~^Mpc {Hq = 100 hkms~^Mpc~^ so that 
a typical value for /i ~ 0.7). 



The auto-correlation in galaxy clustering is used to extract the correlation from a 
galaxy survey, 

«.) = /^(xO^(x.)\ (73) 



P P 

where p and 5p are the homogeneous and fluctuation in the matter denstiy and (...) is the 
average over all points in the sky such that |xi — X2I = s. Figure Isl shows the correlation 



function[3] for 46,748 luminous red galaxies from the Sloan Digital Sky Survey. 

Since galaxy surveys are 3-dimensional, the correlation can be seen in both the ra- 
dial (line-of-sight) distribution, ry, and the angular (transverse) correlations, r±, of the 
galaxies. These are effected differently by the universe's expansion 

ni(^) = (74) 

r±{z) = dA{z)Ae. (75) 
The 3D correlations constrain a combination of these 

Dy{z) = ririi = j^^dAiz?. (76) 

which can then be used to constrain cosmological parameters. 

The Sloan Digital Sky Survey [3] place constraints on Dy at the effective redshift of 
the galaxies in the sample, z = 0.35, 

Ai, . 0.35) . ^^(^ = °-35)v^ ,„) 
^ ' 2c = 0.469 ±0.017 ^ ' 
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Supernova Cosmology Project 
Kowalski, et al., Ap.J. (2008) 




Figure 4: Combined constraints in the Vtm — f^A (left panel) and Vim — w (right panel) 
parameter spaces for the WMAP 5-year CMB data, SDSS BAO and Supernova Cosmology 
Project (SCP) Union supernovae constraints j^. 



The baryon acoustic oscillations, because they have a characteristic size that can be 
determined from an understanding of the physics of recombination, are an additional, 
and complementary, standard ruler to the CMB fluctuations. However while the CMB 
is a standard ruler to z ~ 1100, the BAO provide a standard ruler at the redshift of the 
galaxies used to measure the correlation, at 2; < 1. As such, they provide information, 
complementary to that from supernovae, about the evolution of the expansion history at 
the redshifts when cosmic acceleration is occurring. 

2.4 Complementarity of the distance measures 

In figure |4] we show the complementarity of the supernovae, CMB angular diameter dis- 
tance and BAO datasets, with recent constraints on the fractional energy densities of 
matter and A, and the dark energy equation of statefl]. The three provide consistent, 
overlapping confidence level contours, that in combination indicate a concordance cos- 
mology with ilm = 0.279 ± 0.015, Ha = 0.721 ± 0.015 at 68% confidence level for a flat 
ACDM cosmology, and -0.097 < w + 1 < 0.0142 at the 95% confidence level for a flat 
universe with a constant dark energy equation of state [S] . 

Figure |5] shows the constraints on curvature and the dark energy equation of state if 
the assumption of flatness is relaxed. 
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w w 

Figure 5: Combined constraints in the Q^. ~ w parameter spaces for the WMAP 5-year 
CMB data, SDSS and 2 degree field (2dF) BAO measurements and Supernova Cosmology 
Project (SCP) Union supernovae constraints [5]. 



3 Theoretical avenues to explain cosmic acceleration 
3.1 The Cosmological Constant revisited 



We have already seen in section [T73 that a Cosmological Constant can provide accelerated 
expansion if it dominates the energy density. 

In classical General Relativity, A is a constant of nature able to be added into the 
Einstein-Hilbert action for gravity without destroying covariance 

Seh = I d'x^-^{R - 2A). (78) 

Einstein introduced A to create a static universe (at the time astronomical observations 
had not yet extended out of the galaxy, and so there was no evidence for cosmic expansion) . 
Considering the expansion and acceleration rates today and requiring both to be zero, 

A-l^p = «=0 (80) 

6 a 

gives a constraint relating the curvature and value of A, Rq = c/\/A with k = +1 i.e. a 
closed universe. 

Two quotes from Einstein convey his feelings about the introduction of A 

"/ have again perpetrated something related to the theory of gravitation that 
might endanger me of being committed to a madhouse. " 

"/A is] gravely detrimental to the beauty of the theory. " 
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From a classical perspective, as a constant of nature, it is somewhat philosophical to 
ask why A has the value it does. However from a quantum perspective the value has great 
meaning and its observed value causes a major puzzle. 

If we consider A originating in the vacuum energy, then one can use back-of-the- 
envelope arguments to estimate its expected size. The free-field vacuum energy for a 
mode of frequency uj is huj/2 where uj = y/k"^ + m? with k the wavenumber of the mode 
and m the mass of the particle. The total vacuum energy is then the sum over the vacuum 
energy over all modes and all particles 

particles u \ / j 

where Qi — {—iy^{2j + l) is the degeneracy factor for a particle of spin j, with Qi > 
for bosons and Qi < for fcrmions. p^ac is therefore quartically divergent in k^ax- The 
question now arises, what is kmax^- 

• Pauh in the 1930s calculated the effect of a photon zero point energy on the curvature 

of space and already saw problems. If Xmax is the classical radius of an electron 
~ 10~^^m then the contribution to A = pJJ^^ = a;^„^/87r^ implies Rq ~ Sl/cm and " 
one could not even reach the moon"\ 

• If J?o ~ so that pyac ~ Pcrit '"^ 8 X 10^^*^ Jm~^ ~ lO'^^^GeV^, as we observe, 
then that implies that k^ax ~ lO^^eV (or a cut off scale X^in ~ 10~^m) well below 
scales at which we believe we fully understand physics. 

• If the cut-off is at the Planck scale, then 

Ac^ ^ ^ ~ / ""/"Is ,3 - 10''' Jm-' ^ W'GeV ^ IO^^^ch* (82) 
L^, (1.6 X IQ-^^my 

If one is more conservative about the upper energy limit assuming a cut-off at the 
QCD scale, A ~ lO^^GeV^ ~ lO^^PcHt- Invoking supersymmetry provides a way to 
generate zero vacuum energy while SUSY is unbroken, since bosonic or fermionic 
superpartners exist for each Standard Model particle, with opposite degeneracy 
factor, so that p^ac = Pboson—\Pfermion\ = 0. When SUSY is broken at Tsusy ~ ITeV 
then the contribution to the vacuum energy would be A ~ Pvac^ ^ lO^^GeV^ ~ 

10^^ Pcrit- 

All the options above lead to a phenomenal discrepancy between the theoretical 
quantum value of A and the observed value. It is highly unlikely that a classi- 
cal contribution to A would exactly this quantum contribution precisely enough 
to reconcile theoretical and observed values- this is known as the A fine-tuning 
problem. 

Another puzzle associated with the observed value of A is the Coincidence problem 

- why are the densities of matter and A comparable today? We see that observational con- 
straints suggest that matter comprises 30% of the energy budget today while A makes up 
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most of the remainder. Why are they so similar given that they evolve so very differently 
(A is undiluted by the universe's expansion while matter is diluted as a~^)? 

If A were maybe as little as 10 times larger than the observed value today then it 
would have catastrophic implications for the development of life, because the earlier onset 
of acceleration would suppress or prevent galaxies formation, A < would also have poor 
consequences since the universe would quickly recollapse. In light of this, some authors 
have argued the value of A is influenced by the anthropic principle. Here the value of 
A can vary from one region of space to another, and from one universe to another, and 
the value in our univesre reflects the fact that our universe was able to create galaxies and 
life [SI [71 El [HI Uni [11] . Associated with the anthropic arguments are recent developments 
in the study of string theory vacua. It appears that there could be a very large number of 
possible vacua (> 10^°°) that can arise from the role of gauge field and brane configurations 
in the compactification of extra dimensions [T^ [HI H^l HSl [IZj • 

Raphael Bousso's TASI lectures from 2007 are a good resource [12] for a more in depth 
discussion of the theoretical origins of the cosmological constant . 

3.2 Scalar field dark energy 

In the absence of a robust resolution to the fine-tuning and coincidence problems, alterna- 
tives to A have been considered to give rise to cosmic acceleration. One of the most widely 
investigated ideas is that a new type of matter, a scalar field known as quintessence, is 
driving acceleration. 

The properties of the scalar field can be succintly described by an action, S, 



where d^(f) = d(f)/dx'^ and V{(j)) is a self- interaction potential driving the dynamics of the 
scalar field. 

The energy momentum tensor for the scalar field can be obtained by varying the action 
with respect to the metric 



We will assume that the scalar field is homogeneously distributed in space so that di(j) = 0, 
in which case considering the 00 and ii components of T^^'' give the scalar field's energy 
density and pressure. 



S 






(84) 







1 + 



02 + 2\/ 
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Acceleration is achieved if the scalar field is slowly evolving, 0^ <^ V and dominates the 
energy density. 

The equation of motion for the scalar field can be obtained from the varying the action 
with respect to the field itself, through the Euler-Lagrange equation, 

d(p ^ ddfj^cj) 

which in the FRW metric gives 

+ 3H<i) + V' = (90) 

where V = dV/dcj), and 3if0 is "Hubble drag", an effective frictional term due to the 
expansion acting against evolution in 0. At the terminal velocity under the action of this 
friction term, = and ~ V'/3H so that the condition for acceleration is dependent 
on the form of the potential, 0V^ ~ (V /VY < 1. 

The behavior is effectively a low-energy analogue of the inflaton in the early universe, 
with the advantage that we don't have to address the "graceful exit" problem of how 
inflation is brought to an end, but with the disadvantage that it is far more difficult to 
motivate the existence and dynamical evolution of a scalar field in the low-energy universe 
today. 

Nevertheless, the aim of introducing quintessence is to use £(0) to address the fine- 
tuning and coincidence problems. Let us examine the behavior of for a couple of 
potentials in this regard. 

A power law potential, V = M^"'""/0", gives rise to evolution of the scalar field that 
follows a dynamical attractor, that gives accelerated expansion at late times without 
independently of the initial conditions for the scalar field jl^ . For the scalar field 
energy density to be comparable to the critical density today M ~ IGeV for a = 2 com- 
patible with high energy astrophysics. This would in some way mitigate the fine-tuning 
problem though it still remains to be seen if such potentials are able to be convincingly 
motivated from high energy particle physics. 

An exponential potential, V = rn^e~^'^^^^^ , gives acclerated expansion in the absence 
of matter, but tracks the scahng behavior of matter, i.e. w</, = Wmat otherwise |2Tl|22ll23] . 
This can mitigate the coincidence problem, since the fractional energy density of the scalar 
is a fixed fraction of the background, dominant matter, Q^f, = 3(1 + Wmat)/^^, however 
there then remains the issue of how to generate acceleration at late times. Acceleration 
can be induced through extending the scenario by having a sum of two exponentials, 
including a feature to the exponential potential (to capture and slow the scalar) [21], 
or by adding in a non-minimal coupling between the scalar and matter. The review by 
Copeland et al. [25] provides a comprehensive summary of many possibilities explored. 

Because, for the exponential potential, the scalar energy density has a fixed fractional 
energy density it opens up the possibility that, unlike A, it could make up a significant 
fraction of the energy density at early times (prior to acceleration starting). This is known 
as early dark energy. One can place upper limits on the contribution of at early times 
both at Big Bang Nucleosynthesis and at recombation [26]. At BBN, the scalar behaves 
like relativistic particles with w = 1/3, and boosts the expansion rate as an additional 
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(89) 



relativistic energy density. This alters the expected abundances of deuterium and "^He 
and compared to observed abundances and one finds the constraint flti^{lMeV) < 0.045. 



3.3 General scalar field models 



A variety of extensions to the minimally coupled scalar field model discussed in |3.2| have 
been considered. At its most general we can consider an action of the form 



S = Id Xa/^ 



^7^ + £^(0,X) + £^(0) 



(91) 



where X = g^^dfj_(f)di,(f)/2. The first term represents a non-minimal coupling between (f) 
and gravity, the second is a non-canonical Lagrangian for the scalar field, and the final 
term is a non-minimal coupling between normal matter and 0. We will consider each of 
these possibilities separately. 

3.3.1 Modified gravity 

While A and quintessence address cosmic accleration through a modification to the right 
hand side of Einstein's equation, an alternative might be that GR is modified on scales 
c/H^\ 

f(R) gravity: This is a class of modified gravity models in which the gravitational 
action contains a general function f{R) of the Ricci scalar. 



(92) 



Varying the action with respect to the metric gives the Einstein field equations which now 
contain additional terms on the LHS 

(1 + Jr) R^.u - {R + f) + {g^uO - V^V.) Jr = SnGT^, , (93) 

where = df/dR. 

The Friedmann equation and acceleration equations are modified, 

H' + UHfn = (94) 

--H'fR + a'i + ^T,Hfn + lfR = -^(p + 3P). (95) 
a 6 2 2 6 

The extra terms in the acceleration equation are able to reconcile the observed acceleration 
a > with a universe populated by matter with positive pressure. 

There are a number of conditions on suitable forms of the function f{R)'- i) Jrh > 
so that the dynamical behavior is stable in the high-curvature regime [27] , ii) 1 + Jr > 
for all R so that the effective value of Newton's constant, Ge// = G/{1 + Jr), is positive 
(so that gravitons are not ghost-like), and iii) < so that one recovers GR in the early 
universe. 
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The action (92) is known as the Jordan frame. There exists a complementary, and 



sometimes conceptually simpler, way in which to approach f{R) modifications to GR 
known as the Einstein frame. The two are related by a conformal transformation on 
the metric, 



^Mi. = (1 + fR)g. 



(96) 



so that in the Einstein frame the gravitational action in the usual Einstein Hilbert form of 
GR. The price one pays for this simplification is a non-minimal coupling between matter 
fields and the new metric [2H1 EH], as well as the appearance of a new scalar degree of 
freedom evolving under a potential determined precisely by the original form of the f{R) 
coupling in the Jordan action. 



The gravitational action (92) can be recast into a dynamically equivalent form 
in the Einstein frame 



S 



M2 



+ / d^x^ 



(97) 



where /30 = Mp ln(l + //?), with /3 = a/2/3, and the potential V{(j)) is determined entirely 
by the original form ( 92 ) of the action and is given by 



r(0) 



Ml Rfn - / 
2 (l + fR)' 



{9t 



The Friedmann and acceleration equations in the Einstein frame take the standard forms 
in GR: 



d 
d 



3M2 



p(0) + + v{<P) 



6M2 



p{4>) + 3P(0) + 202 _ 2V{4>) 



(99) 
(100) 



however they contain non-minimal couplings between the scalar and matter, which lead 
to the matter effectively feeling a fifth force. Matter test particles do not move along 
geodesies of the metric g^^, and the density of matter is not purely affected by cosmic 
expansion but also by direct interconversion to 0. This is reflected directly in the scalar 
and matter fluid equations. 



cj) + 3H(j) + V,^ 
P + 3i^ (p + P) 



1 



2M 



2M 



■/3(p-3P) 

I 

■/3(p - 3P) 



(101) 
(102) 



where p = e ^^'f'/^^pp and P = e '^P't'l^vp are the Einstein frame energy density and 
pressure. 
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The Einstein and Jordan frames are wholly equivalent, neither one is more "correct" 
than the other, however one has to think carefully about underlying assumptions made 
the interpretation of observations, for example how redshifts are interpreted, in order to 
match theoretical predictions to observational data. Typically redshifts are interpreted 
assuming minimally coupled matter, i.e. the Jordan frame. 

Extensions of f{R) gravity have also been considered, that involve terms including 
scalar combinations of the Ricci tensor and Riemann tensor, f{R, R^^R^^, R"^^" Rap^^,) 
[22] . An example of this are models in which the action is modified by including the 
Gauss-Bonnet invariant, G = R? — AR^^R^^ + R"^'^" R^p^y, [33J chosen because it does 
not contain any badly behaved "ghost-like" degrees of freedom. However it is difficult to 
find such f{G) models that are consistent with current observational data [M] . 

The DGP model: Another scenario is that gravity is not confined to 4D spacetime, 
and that it is able to pervade extra dimensions. One of the most widely studied extra 
dimension models is the Dvali-Gabadadze-Porratti (DGP) model whose cosmology is 
outlined in [36] . In this theory, matter is confined to a 4D "brane" in a 5D fiat, Minkowski 
"bulk" of infinite volume. The cosmic acceleration we observe is then due to an effective 
4D theory of gravity on the brane. 

The 5D action for DGP is 

M|i?5 , f ,4 / MlR^ 



S = J (fx^^^^ + j rf^xv^^-^ [Ra + CJ^ (103) 

where M„ and i?„ are the metric, Planck mass and Ricci scalar in n dimensions. 
The behavior of gravity in these models transitions at a cross-over scale r^, 

^ (104) 



On smaller distances from the source, r < Tc, gravity appears four- dimensional, with 
$ oc 1/r, while on large scales the graviton's interaction with the fifth dimension becomes 
important and gravity is weakened, with $ cx 1/r^. Though gravity is four-dimensional 
on scales below Tc, it is different from standard General Relativity on distances down to 
the Vainshtein radius, = (r^r^)^/^, where is the Schwarzschild radius for the mass. 
The effective 4D Friedmann equation is 



H^ + ^-e-JH^ + - = —p 105 



where e = +1 is the "self-acceleration" branch of DGP, in which acceleration occurs 
when > rc, and e = — 1 is the "normal branch" in which one requires A to generate 
acceleration. While acceleration can be generated naturally at late times in the self- 
accelerating branch for ~ Hq^, it has been shown that the expansion history does not 
fit observational constraints from SNIa, BAO and the CMB as well as ACDM [23 |3H] • In 
addition, it is believed to suffer from the presence of a ghost-like degree of freedom [39] . 

3.3.2 k-essence 

Quintessence has a canonical kinetic energy term in its Langrangian C = X — V with 
X = —\g'^''dn4>d^(j), and relies on a slow evolution of the scalar field down a fiat potential to 
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drive late-time acceleration. An alternative scenario is that accleration is brought about 
because of a modification to the form of the kinetic Lagrangian, so-called k-essence 

HH 1321 US] • This concept is a low-energy analogue of "k-infiation" to explain inflation 
in the early universe |33]- 

This generalized scalar field can be described by an action of the form 



S = J (fx^^C{X,(P) (106) 

which can contain terms including non-trivial functions of X and combinations of X and 
(f). A possible motivation for actions of this form comes from considering low-energy 
effective actions from string theory in which higher-order derivative terms are important 
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S = j d'x^g [b,{<P)R + BfmWct^f - a' [c^^^ B^^\ct^){Wct^Y + ...] + 0{a'Y] (107) 

where is the dilaton field related to the string coupling by Qs = g"^-, is related to the 
string length scale, A^, a' = Xs/2tt. In the weak coupling regime Qs <^ 1, gives couplings 
Bg !=s B^^ ^ B'^^ ^ e~'^. On doing a conformal transformation g^i, = Bg{(j))g^i, the model 
has an Einstein frame action that contains a k-essence style scalar Lagrangian 



S = I di^Xy/^ 



Ml 



2 



3 fdlnBA 5, 



^R + K{(t))X + L{(I))X^ + ... 

(0) 



(108) 



m = 2c«a'i?i^)(0). (110) 
With a scalar field redefinition this transforms into a Lagrangian of the form 

= /(0)(-X + X2) (111) 

for which the equation of state is 

1 -X , , 

= (112) 

and acceleration is obtained for X < 2/3. 

In quintessence the minimum equation of state the scalar field can have is w = —1 
when X = 0, in contrast k-essence models can exhibit phantom like behavior with w < —1. 
Consider a general Lagrangian of the form C = K{X) — V{(j)), by varying the Lagrangian 
with respect to the metric the energy density and equation of state are given by 

p = 2XKx-K + V (113) 

K-V 2XKx 

w = = — 1 H 114 

2XKx - K + V 2XKx - K + V ^ ^ 

where Kx = dK/dX. For Kx < phantom behavior arises. 
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3.3.3 Non-minimally coupled dark energy 



A theory in which the scalar field is non-minimally coupled to matter can be described 
by the following action 



+ 5^5j[e2"^(^V,vI/^], (115) 



where g^i, is the Einstein frame metric, is a scalar field which acts as dark energy, and 
\E'j are the matter fields. The functions aj{(f)) are coupling functions that determine the 
strength of the coupling of the jth matter sector to the scalar field. 

The Friedmann equation and fiuid equations for the scalar field and matter are 



1 V 



a, (4,) 



Ph 



4> + 3H<j) + V'i(l)) 
Pj + 3(l + ti;j)ifpj 



5^a;(0)(l-3«;j)e 



Ph 



(116) 
(117) 
(118) 



where a dot represents a derivative with respect to the Einstein time coordinate, t, and 
primes are with respect to 0, and H = a/a is the Einstein frame Hubble parameter. 
We can solve equations (116)- (118) in the scenario in which CDM alone is coupled to 

28] dimensionless quantities: 



the scalar, by studying Copeland et. al.'s 



X 



V6HM„ 



y 



A 



MpV 
V 



VV" 



:il9) 



We also define a dimensionless variable to describe the coupling to matter, 

«'(0)Mp 



/3 



(120) 



Rewriting in terms of the dependent variable = ln(a) yields 



dx 
dN 

dy 
dN 
dX 

dN 



V 6 , n 3 , n 

3x + —Xy + -x{l + x^ 



+Cil-x'-y') 
— ^^xy + -y{l 
-VQX\r-l)x. 



2 2\ 

X -y ) 



(121) 
(122) 

(123) 



In these equations, F, C and a are understood to be the functions of A obtained by 
eliminating in Eqs. (119) and (120). After the equations have been solved, the matter 
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density p = Pc + Pb (neglecting radiation in the matter dominated era) can be obtained 
from the equation 

^ rPe" = 1 . (124) 



2 I 2 

X + y 



Note that the effective total equation of state parameter Wes, defined by a{t) oc t2/3(i+Woff)^ 
is simply given by 

(125) 



from Eq. ( |124[ ). 

The fixed points of this system are the solutions of the equations dx/dN = dy/dN = 
dX/dN = 0. 

A specific example is provided by the exponential potential V{(f)) = Vqc'^'^/^'^p where A 
and Vq are constants, and by C(0) = C, a constant. The corresponding coupling function 
a{(j)) is then linear: 

a(0) = -l3C<p/Mp. (126) 



The fixed point followed in the matter dominated era is described by 

^2C 



{x,y) 



3 

~9~' 



,0 



(127) 
(128) 



It doesn't give rise to acceleration but deviates from CDM evolution (with Weff = 0) 
depending on the coupling strength. 

A second attractor is independent of the non-minimal coupling and gives acceleration 
for < 2 



ix,y) 



A 

\/6 



A^ 



2\ 1/2' 



6 



(129) 

WeS = • (130) 

For a power law potential, V oc the CDM dominated era has the same attractor 



3 ■ 



(127) while the accelerated era attractor is different 



ix,y) 

^eff 



[0,1] 
-1 



(131) 
(132) 



independent of the coupling, C, and the precise power law form of the potential, n. 

Constraints on the non-minimal coupling strength, C in light of recent CMB, BAO, 
supernovae and galaxy distribution data are \C\ < 0.13 for an exponential potential and 
—0.055 < C < 0.066 for the power law potential at the 95% confidence level 
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4 Using the growth of structure to understand the 
origin of cosmic acceleration 

As discussed in section |2j cosmo logical constraints on the properties of cosmic acceleration 
have, to date, primarily come from the homogeneous background expansion, and asso- 
ciated geometrical distance measures of the CMB angular diameter distance, the BAO 
acoustic scale and supernovae luminosity distances. However a projected increase in the 
number, depth and breadth of large scale structure surveys offer the opportunity to also 
constrain the origin of acceleration by investigating its affect on the growth of cosmic 
structure (the growth of over-densities and peculiar velocities that seed galaxies and clus- 
ters of galaxies). By contrasting expansion history and growth of structure measurements 
one can investigate evidence for dark energy perturbations and modifications to gravity 

imilHlliniEniEIlESlESlIM!. 

4.1 Linear perturbation theory 

In order to study the growth of inhomogeneities we need to consider a perturbed, inho- 
mogeneous metric and energy momentum tensor. A general perturbed FRW metric is 
given by 

ds'^ = -(1 + 2AY)dt^ - aBYidx'dt + a^{5i] + 2HLY6ij + 2HTYij)dx'dx^ (133) 

where Y,Yj and Yij are unit scalar, vector and tensor perturbations respectively, and A,B, 
Hl and Ht are arbitrary time and space dependent scalar perturbations. 

Though it is certainly possible to perform the analysis in this general "gauge ready" 
approach it is somewhat time consuming. Instead we will use an intuitive gauge (choice of 
A,B,Hl and Ht) known as the conformal Newtonian gauge in which A = , B = Ht = 0, 
Hl = $. In this gauge $ reduces to the standard Newtonian potential (that obeys 
Poisson's equation) on scales much smaller than the horizon, and for GR, in the matter 
dominated and accelerated eras, $ = — \E'. 

As well as perturbing the gravitational metric, we also need to describe the perturbed 
RHS of Einstein's equation, the perturbed energy-momentum tensor 

7^0 = -(pit) + 5pix,t)) = -Pirn + Six, t)) (134) 
TO = (pit) + P{t))vY, (135) 

T; = P{t)6'^+6P{x,t)6'^Y + ^{p{t) + P{t))a{x,y)Y; (136) 

where 6 is the fractional energy density, v is the peculiar velocity, and 6P and a are the 
isotropic pressure and anisotropic shear stress perturbations respectively. The isotropic 
shear stress is sometimes rewritten in terms of the speed of sound of a wave propogating 
in that matter species 

6P = cl6p. (137) 
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The perturbed Einstein equations relate the metric and matter peturbations 



= 47rGa2^p,Ai (138) 

i 

+ = -l2TiGa^^{pi + Pi)ai (139) 

i 

The first equation is the Poisson equation with Aj = (5j + 3if(l + Wj)^) the density per- 
turbation in the matter's rest frame. The second reflects that in regimes where relativistic 
species (and hence anisotropic shear stresses) are negligible then $ = — 

The comobing Hubble length, c/aH, gives a measure of the region that is able to 
communicate at a given time. Comoving modes for which k > aH are "subhorizon" while 
k < aH are "superhorizon" . 

Energy momentum conservation, T'^'^-^j, = 0, for the inhomogeneous matter gives us 
the perturbed fluid equations. In terms of conformal time, 6 = dS/dr, and H = dlna/dr, 
for cold dark matter 

5c = -ikvc-3^ (140) 
= -q-iv^-ikm. (141) 

Combined these give a second order differential equation in 5, 

6c - H5c + + 3$ + A;^^ = 0. (142) 

On subhorizon scales the time derivatives of the potentials are subdominant to the term 
in k"^ 

5c + Wc - AttGo^ Y1 Pi^i = 0- (143) 

i 

During accelerated expansion the second term, a Hubble damping term, drives suppression 
of growth in 6c- 

For general matter 



-{l + w){ikv + 3^) -3'H{cl-w)6 (144) 

^2 



which gives 



-n(l-3w)v V —ik6-ikm + ika. (145) 

^ ' l+w l+w ^ ' 



S + nS- {iTTGa^ piSi - clk'^S) = 0. (146) 



If > it reduces the effective strength of gravity on scales k < oHjcs and suppresses 
the growth of the overdensity. For < catastrophic growth of the inhomogeneity can 
occur. This instability has been particularly discussed in the context of non-minimally 
coupled dark energy-dark matter models [55l 
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One might naively think that = = P/p, where is the adiabatic sound speed 



This is the case for perfect fluids such as CDM, baryons and photons. For scalar field dark 
energy, however, this would imply that a negative equation of state would often lead to 
catastrophic collapse of dark energy over-densities, and the formation of dense 'nuggets' 
of dark energy on subhorizon scales, something we do not observe. 

In general a fluid need not be perfect in this way, and the pressure perturbations 
are not wholly specified by the background expansion history. This leads to a disparity 
between and because of "entropy perturbations, 

5Pen ^ {cl-cl)5. (148) 

For minimally coupled quintessence one finds = 1 irrespective of its equation of state. 
For a more general, k-essence model depends on the form of the kinetic action 

P,X + '^^^,XX 

which opens up the possibility of < and an associated instability leading to uncon- 
trolled growth on scales below Cg/aH. 

In modified gravity models, the Einstein equations are modified while the fluid equa- 
tions remain unchanged. A common parameterization introduces two scale and time- 
dependent functions, Geff{k,a) and R{k,a), 

e<^> = 47rGe//a'5^Mi (150) 

i 

^ + = _i27rG'e//a2^p,(l + w;)^. (151) 

i 

The first describes a modified Poisson equation in which the gravitational potential re- 
sponds differently to the presence of matter, with an effective Newton's constant, Ge//, 
while the second allows an inequality between the two gravitational potentials, even at 
late times when anisotropic shear stresses are negligible. In the DGP model [22] discussed 



in 3.3.1 



R = V (153) 



1-^ 



where 



/3 = 1- 2rcH I 1 + I . (154) 



3iJ2 
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For f{R) theories, Ge// = + fn) while i? is a more involved function obtained by 
integrating the Einstein and fluid equations [571 EH] • 

By measuring the evolution of the gravitational potentials, and their relationships to 
the growth of overdensities and peculiar velocities we can search for signatures of dark 
energy innhomogeneities, or modifications to gravity associated with cosmic acceleration. 



4.2 Observations of large scale structure 

For a range of different observations, for example the Integrated Sachs Wolfe (ISW) effect 
in the CMB, galaxy number counts and weak gravitational lensing we discuss below, the 
measurements are based on the two-point correlation function between fluctuations. 

Under the Limber approximation, the 2D angular power spectrum for the correlation 
between two fields, X and Y, is 

/•Xoo 

Cr = / -4Wx{x)Wy{x)Tx{h,x)Ty{h,x)AUki), (155) 
Jo X 

where Wx is the window function which gives the redshift window for the field X, and Tx 
is the transfer function that relates the inhomogeneity in X today to its primordial value, 
and A'^{k) is the dimensionless primordial spectrum of curvature fluctuations, ki ^ l/x 
where I is the multipole moment (analogue of Fourier mode on the spherical sky) and x 
is the comoving distance. 

The statistical uncertainty in the power spectrum is given by 



(2/ + l)fsky 



Ci + 



n. 



eff. 



(156) 



where fsky is the fraction of the sky covered by a survey, cr^ is the variance for a single 
measurement, and Uefj is the effective number density per steradian of measurements 
used in the analysis. The first term in the square bracket represents "cosmic variance", 
the intrinsic sample variance coming from only having a finite number of independent sky 
samples to measure correlations. This term dominates on large scales where the number 
of independent samples is small. The second term is the statistical "shot noise" error 
coming from intrinsic variations in the objects being measured and instrumental noise. 



4.2.1 Integrated Sachs Wolfe (ISW) effect 

In section [2] we discussed how the angular diameter distance to the last scattering surface 
of the CMB provides an important distance measure to constrain dark energy properties. 
The large scale CMB temperature correlations provide additional information that can 
also be used to constrain dark energy however. On scales below the sound horizon at last 
scattering, the physics of Thomson scattering and electron/ baryon interaction dominate 
the CMB temperature correlations leading to the acoustic peaks. On scales larger than the 
sound horizon only gravity and the cosmic expansion history determine the correlations. 
Of particular interest in the context of dark energy is how the recent onset of cosmic 
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acceleration affects the correlations in the CMB on the very largest scales observable, the 
Integrated Sachs Wolfe (ISW) effect. 

The ISW effect arises because of the effect of a time-changing gravitational potential on 
the CMB photon energy. As the photons stream to us today they traverse gravitational 
potential wells, gaining energy as the photons fall into the well, and losing it as they 



climb out. In the matter dominated era, solving (143) one finds 5c oc a, and hence from 



(150) that $ = is constant. This means that there would be no net change in the 
photon energy after moving through the constant gravitational well. During the period 
of accelerated expansion however the suppression of growth in the overdensities causes 
gravitational potential wells to decay so that photons gain more energy on their descent 
into the well than they lose on their climb out. This gain in energy is translated into a 
boost in large scale correlation amplitude in the CMB power spectrum. 

The CMB correlation comes from integrating along the photon's history since last 



scattering so that all redshifts are sampled, in (155) therefore Wisw = 1 and the transfer 



function is dependent on the rate of change of the gravitational potentials 

Tjsw = e""'-"°" f-l> + ^) , (157) 



where X is the transfer function of X, normalized so that X'^{k,x) = X'^{k,x)^Tzi^)- 
Treion IS the reiouizatiou optical depth, taking into account that some CMB photons are 
scattered by free electrons created when the universe was reionized by the first stars 
z ^ 10. 

The ISW signal, being a large scale effect, is dominated by cosmic variance. It has 
been measured to cosmic variance limits by the WMAP survey 



4.2.2 Weak gravitational lensing 

When light travels through a gravitational potential well its path is distorted much as 
it would be through an optical lens. There are two primary effects of the lensing: the 
image can be magnified and the image can be distorted, or undergo "shear", which can 
be, in a simple way, thought of as a circular source having an elliptical image. The lensing 
can be thought of as a remapping of the 2D angular distribution of the source's surface 
brightness, f^, to the observed image, f"^"^, through a "distortion matrix", Aij 

r'\e,) = riA.e,) (158) 

with 

'^'^ = df,-[ -72 l-^ + liJ- ^^^^^ 

Here k is the convergence which is related to the magnification, and 7 = 71 + i'~f2 is the 
shear. 

The convergence is given by 

K = k\<^ + ^)W,{x.Xs)dx (160) 
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where Xs is the comoving distance to the source. The window and transfer function in 
(155) are 

= -y($ + vl/) (161) 

W^{x.Xs) = xTdxHix')^^ (162) 
Jo X 

where n{x) is the number density distribution of lensing galaxies along the line of sight 
normalized such that J^" n{x) = 1- 

Lensing is sensitive to both dark energy's effect on cosmic geometry, through the 
expansion history's effect on the window function, and the growth of structure, sensitive to 
both the expansion history and the relation between gravity and matter inhomogeneities, 
through the transfer function. 

Lensing surveys such as those from the Canada France Hawaii Legacy Telescope 
(CFHTLS) ISni Ell E2] and the HST COSMOS survey [Ml El] have provided the first 
applications of lensing data to test the nature of dark energy. Planned surveys such as 
the Dark Energy Survey (DES), a ~5000 square degree survey with galaxy survey density 
of ~ 15 per square arcmin and mean redshift ~ 0.7, and the Large Synoptic Survey Tele- 
scope (LSST), which will cover ~ 20,000 square degrees of sky, with ~ 30 galaxies per 
square arcmin and mean redshift 1.2, hope to provide richer and higher precision lensing 
data with which to test the origins of cosmic acceleration. 

Systematic errors in weak lensing measurements include intrinsic correlations (intrin- 
sic alignments) in galaxy shapes, uncertainties/ errors in photometric redshift estimates 
of galaxies and errors in shape measurement arising from anistropy of the point spread 
function (PSF) of the detector, atmospheric distortions and pointing errors amongst oth- 
ers. 



4.2.3 Cross-correlation of large scale structure observations 

The ISW and lensing observations directly probe the properties of the gravitational po- 
tential. By contrast, measurements of the distribution of luminous matter provide only a 
proxy for the gravitational potential. To recover the distribution of all matter, both dark 
and luminous, one needs information about how the distributions of dark matter and lumi- 
nous matter are related, through a potentially scale and redshift dependent bias. However 
by using cross-correlation of galaxy, ISW and lensing variables one can reduce the sensi- 
tivity to bias With regards to understanding cosmic acceleration it has been shown 
that cross-correlations could provide a powerful tool with which to distinguish between 
modified gravtiy and dark energy by measuring Geff/G and the relationship between the 
two Newtonian potentials ^/^ [HSlESlEIlEHlEaiiniinilZSllISllTilllS]. 



5 Final remarks 

These lectures have discussed how gravity on cosmic scales relates the evolution of space- 
time scale and curvature to the properties of matter contained within it. Current observa- 
tions of geometric distances to Type la supernovae, baryon acoustic oscillations and the 
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Cosmic Microwave Background last scattering surface have provided consistent, strong 
evidence that the universe is undergoing accelerated expansion. This is seemingly at dis- 
cord with the a universe in which gravity is described by GR and matter has positive 
pressure. We have considered a spectrum of possible theoretical origins for this cosmic 
acceleration: a perfectly homogeneous energy source with negative pressure, a "cosmolog- 
ical constant" ; a near-homogeneous, potentially dynamically evolving new form of matter, 
" dark energy"; or that gravity deviates from GR on cosmic scales, "modified gravity". 
We have discussed how the combination, and comparison, of observations of both geo- 
metric distances and the growth of large scale structure, such as the ISW effect, weak 
gravitational lensing and galaxy number counts, could potentially give us an insight into 
which of the theoretical possibilities is the origin of cosmic acceleration. Upcoming sur- 
veys in next few years should provide excellent data with which to address this central 
issue facing cosmology and particle physics today. 

6 Further reading 

The following reviews are good resources for further reading on dark energy and cosmic 
acceleration: 

• "Dynamics of dark energy" [25], Copeland, Sami and Tsujikawa, 2006 

• "Dark energy and the accelerating universe" [76] , Frieman, Turner and Huterer, 2008 

• "Approaches to understanding cosmic acceleration" [77], Silvestri and Trodden, 2009 
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